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A class of two-species [three- states) bimolecular diffusion- limited models of classical particles 
with hard-core reacting and diffusing in a hypercubic lattice of arbitrary dimension is investigated. 
The manifolds on which the equations of motion of the correlation functions close, are determined 
explicitly. This property allows to solve for the density and the two-point (two-time) correlation 
functions in arbitrary dimension for both, a translation invariant class and another one where 
translation invariance is broken. Systems with correlated as well as uncorrelated, yet random initial 
states can also be treated exactly by this approach. We discuss the asymptotic behavior of density 
and correlation functions in the various cases. The dynamics studied is very rich. 
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I. INTRODUCTION 



Non-equilibrium statistical mechanics has wittnessed recently a resurgence of interest. Though over 50 years old, 
the field is still in its infancy. Powerful concepts and tools are being developped, and yet much progress remains to be 
done. The understanding of classical stochastic many-body systems is of relevance to a wide class of phenomena in 
physics and beyond. In this context, a class of models describing diffusion-limited reactions plays an important role 
[Q. The natural language to describe the stochastic dynamics of N classical bodies is that of the Master equation. 
Formally, the dynamics can be coded in an imaginary time Schrodinger equation, where the Markov generator plays 
the role of the Hamiltonian (see e.g. |^J^ and references therein). In the past, various representations in terms of 
spins, fermions or bosons have been used depending on the physics being emphasized. A powerful method that 
relies on bosonic field theory and the renomalization group has been applied by Cardy and collaborators to deal 
with low density systems in arbitrary dimensions (see e.g. |^). In one spatial dimension alternative approaches have 
been proposed which take into account the hard core of the classical particles (see e.g. and references therein). 
Since the pioneering work of Glauber on the stochastic Ising model y] various generalizations and extensions have 
appeared (e.g. |^,^) . A general approach has been proposed by Schiitz Q. The latter investigates the most general 
class of single-species models of bimolecular diffusion-limited reactions that can be solved exactly. Upon imposing 
constraints on the available manifold, the equations of motion of all correlation functions close and, in that sense, the 
dynamics is completely soluble. Via a duality transformation, Schiitz further shows that on the 10-parametric soluble 
manifold, the spectrum of the stochastic Hamiltonian coincides with that of the XXZ— Heisenberg model (see also 
[H). In Reference Fujii and Wadati extend Schiitz's ideas to the s— species models of bimolecular diffusion-limited 
reaction processes. These authors derive the general constraints that allow for the equations of motion of correlation 
functions to close and, similarly to the single species case, introduce a dual Hamiltonian with identical spectrum. 
They further note that in the general multi-species case, the constraints of solubility (in the sense given above) do 
not seem to imply a simple relationship to integrable quantum Hamiltonians. 

The general s— species bimolecular reaction-diffusion processes are characterized by (s -I- 1)^ — (s -f 1)'^ independent 
parameters (reaction-rates) and we have to impose 2s^ constraints to close the hierarchy of the equations of motion 
of correlation functions. As few exact and explicit results for the dynamics of multi-species processes are available (in 
particular in dimensions c? > 1, see the discussion at the end of section HI) , we decided to investigate in some details 
and generality the soluble two-species bimolecular diffusion-limited reaction systems. In this paper, we focus on the 
two-species problem (s = 2) and obtain, in arbitrary dimension, exact results. A particular physical application of 
this work to a three-states growth model will be presented elsewhere [ p^ 

The paper is organized as follows: the remaining of this section will be devoted to definitions and notations. In section 
II, the equations of motion of the density and two-point correlation functions are derived. The constraints that ensure 
the solubility of the problem are explicitly identified. We classify the soluble manifolds which will be investigated in 
the sequel. In the first part of section III, we study the Fourier-Laplace transform of the density in the soluble case 
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(on a 56— parameters manifold). In the second part, we compute on two manifolds the exact expression of the density 
in arbitrary dimensions. We provide the asymptotic behavior of the latter for three different initial conditions. At the 
end of section III we discuss the relationship between our results and the solution of some models solved exactly in 
dimension d > 1 fll|-[l6|] . In the first part of section IV, we give the exact dynamic form factor for an homogeneous and 
uncorrelated (yet random) initial state. In the second part, we compute, in arbitrary dimension, the exact two-time 
two-point correlation functions for random (homogeneous) uncorrelated as well as correlated initial states (we discuss 
the sensitiveness of the system to the presence of initial correlations) . Section V is devoted to the study of the 
instantaneous two-points correlation functions on a manifold of translationally-invariant models. We first deal with 
the one-dimensional case, which is investigated for random uncorrelated as well as correlated initial states. Further, 
we consider the higher dimensional case with random (yet homogeneous) and uncorrelated initial states. The last 
section, is devoted to the conclusion. 

For clarity and brevity's sake, some definitions, as well as some technical details, are given in the appendices. 

Consider an hypercubic lattice of dimension d with periodic boundary conditions and N sites {N = i''), where 
L represents the linear dimension of the hypercube. Further assume that local bimolecular reactions between species 
A and B take place. Each site is either empty (denoted by the symbol 0) or occupied at most by one particle of type 
A (respectively B) denoted in the following by the index 1 (respectively 2). The dynamics is parametrized by the 
transition rates r^^^*, where (3i , (32 , , ^4 = 0, 1, 2: 

V(/3i, /52) ^ if33,(34) , rff^t ■■ Pl+P2-^P3 + (34 (1) 

Probability conservation implies 

(/3i,/32)/(/3;,/3^) 



with 



r^^^:>0, V(/3i,/32)^(/33,/34) (3) 



For example the rate Fjf corresponds to the process A + A — > A + B, while conservation of probability leads to 

pll _ /ylO 1 pOl 1 pOO 1 p02 I p20 i p21 _i_ pl2 , p22\ 

1 11 — —[^ 11 -I- i 11 -I- i 11 -I- i 11 -I- i 11 + 11 + 11 + i 11 j 

The state of the system is determined by specifying the probability for the occurence of configuration {n} at time 
t. It is represented by the ket 

\Pit))=J2Pi{n},t)\{n}) (4) 
{»} 

where the sum runs over the 3^ configurations {N — L"*). At site i the local state is denoted by the ket \ni) = (10 0)'^ 
if the site i is empty, \ni) = (0 1 0)"^ if the site i is occupied by a particle of type A (1) and jn^) = (0 1)"^ otherwise. 
It is by now well known that the Master equation governing the dynamics of the systems can be rewritten as an 
imaginary-time Schrodinger equation: 

||P(t))=-ff|P(t)), (5) 

where H is the Markov generator, also called stochastic Hamiltonian, which in general is neither hermitian nor 
normal. The construction of the stochastic Hamiltonian H from the master equation follows a known procedure (see, 
e.g. i|,i0i]). We define 11,101] the left vacuum {x\ : 

(xl^E(MI (6) 
{»} 

Probability conservation yields the local equation (stochasticity of H) 

{X\H = E(x|i/™,™+e» = =^ {x\Hm,m+e^ = 0, (7) 

a=l,...,d m 

where e", 1 < a < d, designates, in cartesian coordinates, the unit vector along the a-direction. 
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The two-species local Markov generator acts on two adjacent sites, i.e., 
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where the same notations as in reference Q have been used. Probability conservation implies that each column in 
the above representation sums up to zero. Locally, the left vacuum {x\ has the representation 



(xl (1 1 1) «) (111) = (11 11 11 111). 



(9) 



The action of any operator on the left-vacuum has then a simple summation interpretation. This observation will be 
crucial in the following computations. Below we shall assume an initial state |i^(0)) and investigate the expectation 
value of an operator O (observables such as density, etc.) 



{0){t)^{x\Oe-"'\Pm 



(10) 



II. DYNAMICAL EQUATIONS OF MOTION 



Exploiting the properties of the left vacuum {x\ and denoting by rif,/3 G {0,1,2} the occupation of site i by a 
particle of type P, we derive below the equations of motion of the density and two-point correlation function. For 
/3 = 0, n° denotes the empty state at site «, i.e.. 



(11) 



We evaluate the action of H on the operators nf, nf , taking into account the local nature of the Markov generator: 



(12) 



As an example, the first term in the above yields 
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where the use of equation (g) is required to substitute for F}}, F^Jq, rj^. 
connects the one-body initial operator to a two-body expression. 
The equation of motion for the density becomes (at site m) 



(13) 

As expected, the stochastic Hamiltonian 
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In order to determine the second moments, we also need to evaluate the following terms 



For the sake of illustration, the first term above yields 
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Notice that the evolution operator connects a two-body operator to a two-body expression. 

To compute the two-point correlation functions, we have to distinguish the sites that are nearest neighbors from 
those that are not. If the sites m andn arc not nearest neighbors {dist{m, n) > 1), the equation of motion reads 



d_ 
dt 
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while if the sites are nearest neighbors, we have 
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The equations of motion of n— points correlation functions are obtained in a similar way, with help of ( p^ , |T^) and 
( p^ , . As is well known, the equations of motion of classical (or quantum) correlation functions constitute an open 
hierarchy which is not soluble in general. However, if we impose on the 3** — S'^ — 72 bimolecular transition rates 



involving two adjacent sites, the following 16 contraints (see appendix A for the definitions (141,142)) 
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the equations of motion of the density and two-point correlation functions (and all multi-points correlation functions) 
become closed. It is worth emphasizing that when the hierarchy closes at the lowest level, i.e., at the level of the 
density, the equations of motion of all higher correlation functions also close. This is a remarkable property. 
In the expressions (|l^) and (^9|), the rates F"^ have not been made explicit for brevity. 

A general diffusion-limited two-species reaction model is defined on the manifold, Vpar — — {r"^}|a,/3 € 

(0,1, 2)|, which has here 3^ — 9 — 72 independent parameters. Let us denote by Vsoi the restriction of Vpar on 

the (72 — 16 = 56 parameters) manifold defined by the additional constraints (p^: Vsoi = Vpar H (p^). The latter 
represents the manifold on which the equations of motion of the correlation functions are closed, i.e., the soluble 
manifold. We can further require translation invariance, i.e., {nl^n-'^^^^^^^){t) ~ (nQn|'^|)(t) = Gi^^\{t), ^r^t G {A,B)) 
and in particular {n^n^){t) = {n^n^)(t) = G^^^,^t). Imposing the above conditions in equations ( p7| ) and(p^ and 
taking into account the conditions of solubility (19), we arrive at the manifold Vtr an si -invar, the restriction of Vsoi 
on the translation invariant soluble dynamics. Notice that Vtransi~-invar{d) — Vsoi H V'{d), where V'{d) ~ ^E^'' = 
^ba. pab _^ pab + ^[jd + (d - 1) = i^^" + F|" + Cl^d + A'^{d - I); + Ff + Al{d~l) + C^d^ F^'' + F^" + (d - 
1) + Agd; Hf + + (Cf + B^d + {A1 + D\){d - 1) = Hf" + H^'' + (A? + D^)d + {B\ + )(d - 1); A'^ + D\ = 
B^ + Cf; G\'' + C^d + A\{d-1) ^ Gf + C^{d - 1) + A\d; Gf + B1{d - 1) + D^d ^ Gl'' + Bfd + D^{d - 1); B^ = 

-D"; A2 — C^j. Therefore this manifold has 72 — 16 — 9 = 47 independent parameters. In practice, however further 
constraints may be required for the computations to be accessible. With this remark in mind, we define the manifolds 

V" = + C| = = = Bf + = B^ = = 0^ and V" = = Bf^ = G^, = D"^ = G^ = G\ = Gf = 

Gf = H-^b = 0\n = l,2y 

Summarizing the cases that we will discuss in this paper: 

i) For the case where translation invariance is broken, we shall compute the exact density on the manifold Vi , which 
has 72 — 16 — 6 — 50 independent parameters, 

Vi = nVsoi n V" (20) 

ii) For the translation invariant case, we shall evaluate both, the density and two-point correlation function exactly 
on the manifold V2 (d) , 

V2{d) = Vol n V'{d) n V" = Vtransl-^nvar{d) H V" (21) 

which has 47 — 16 = 31 independent parameters. 

To conclude this section, it is worth noting that there are few cases in which the open hierarchy of equations of 
motion can be solved analytically. The class of single-species one-dimensional models for which the evolution operator 
can be cast into a free fermionic form is an important example. However, the procedure of free " fermionization" 
cannot be applied to higher dimensions and/or multispecies problems, contrary to the method followed here. 



III. THE DENSITY: GENERAL DISCUSSION 



In the first part of this section we compute exactly the Fourier-Laplace transform of the density on the 
56— parameters manifold Vsoi, which is, as are correlation functions, directly related to light scattering measurements 
in real reaction-diffusion systems ||l^-|l^ . The computation on the most general soluble manifold is here manageable 
because the linear differential difference equations governing the dynamics give rise to a 2 x 2 matrix, the properties of 
which can be studied analytically. For higher order correlation functions and/or for the s-species case, with s > 2, the 
problem is however technically much harder (we shall come back to the general case in a future work). In the second 
part of this section we provide the density of species A and B in space and time, both in the translation invariant 
case and in a situation where translation invariance is broken. On the manifold (Vpar H Vsoi) 3 {Vtransi-invar H Vsoi), 
we have 

J^int) = (^S + C^,)d+{ni){t){A'i + G^2)d + Y.{^2{^t+e-)i^) + (n;^_,.) (t)) 

a 

+ (0(i)(i?? +i?2")rf + E +^?("™-e">W) (22) 
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= (4 + C'o)d + {nf^){t){B\ + Dl)d + ^ (i?2'(n^+e=) W + (0) 



Let us first consider tlie most general soluble case which is characterized by the set of equations ( |22| , p3[ ) . The solution 
of ( p^ , p3| ) is split into the solution of the homogeneous system {n^)h{t) {{n^i)h{t)) and a function fA{t) that 
takes into account the inhomogeneity, i.e., 

jjA{t) = (Ag + CS)d + fAmAI +A-^+ CI + C^)d + fB{t){Bl + B-^ +DI + D^)d, (24) 

^/B(t) = [A'o + C'o)d + fsmsl + Bl + D\ + Dl)d + fA{t){A\ + Al + Cl + Cl)d, (25) 
We introduce the Fourier transforms of {n^-^^)h{t), i.e., 

pel.B.Z. in 

where the sum on p runs over the first Brillouin zone (I.B.Z.). The solution of the homogeneous problem in Fourier 
space reads 

where A4ij{p), G (1, 2) is a s x s = 2 x 2 matrix with the entries 

Mi,i{p) = [Al + C^)d + E (A^e*P-^° + C1e-'P-^-) 

a 

Mi,2{p) = [Bl + D^)d + {B^e'P"" + D'^e-'P-"-) 

a 

M2,i{p) = + Cl)d + Y {A\e'P-'- + C\e-'P-'-) 

M2,2{p) = (Sj + D'2)d + Y {Ble'P"- + D\e''P"-) (28) 

a 

The eigenvalues of the matrix Ai , which represent the inverse relaxation times of the system, control the asymptotic 
behavior of the density. 



7Wi,i(p) + M2,2{p) , y/{Ml,lip) - M2,2{p)f + 4A^i,2(p)A^2,i(p) 
A±(W = ^ ± ^ (29) 

It has been shown in considering the one-dimensional alternating-bonds Ising model obeying Glauber's dynamics pc| ] 
that the relaxational eigenvalues of the analog of the matrix M. allow to identify the critical (but non-universal) 
behavior: it is determined by the long wave-length p modes of the analog of the acoustic A_ branch. One-dimensional 
alternating-bonds (Ji > J2 > 0) Ising model, with Glauber's dynamics, exhibits a non-universal critical dynamical 
exponent z = 1 + ^ 

In the sequel we shall need the zero- momentum 2x2 matrix M.{p = {)) = M. (0) 

7Wi,i(0) = {Al + AI + Cl + Cl)d 
A^i,2(0) = {Bl + B'i + DI + D^)d 
M2a(0) = {A'i + + + C!i)d 
M2AO) = {B\ + Bl + D\ + Dl)d 

whose eigenvalues we shall denote, for short, 
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7± = A±(P = 0) 



Notice that at p = 0, TrM{0) < and detM{Q) > 

We are now in a position to compute the Fourier-Laplace transform Sq'^{p,u)) of the density : Sq'^{p,co) = 

We consider initial states (n^-^)(0) = Er„(«m'^)(0)e-*P"'-"(^ («m'^)(0) = EpeiB.z.(^p '■^)(0)e'^-™) and assume 
that the matrix A^(p) is regular (i.e. detAi{p) ^ 0). We shall distinguish four cases. 

1) We assume that the matrix M.(j)) is diagonalizable, i.e., A_)_(p) ^ A_(p), but not triangular, M.\^i(j>) ^ 0, and 
obtain 
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(30) 



«(0)) 



(31) 



Notice that the inhomogeneous part of the equations of motion give rise to a zero-momentum contribution which we 
shall omit hereafter. As expected the poles in the w-plane occur at the relaxational eigenvalues. 

2) Next consider the case where the matrix M. (p) is non-diagonalizable and non-triangular: X{p) = A+ (p) = A_ (p) = 
MiAp)+M2.2{p) and 7Wi,i(p) ^ M2,2{p) ■ We can compute e^^^)* with the help of a Jordan decomposition of the 
matrix M.{p), namely, 



(32) 



where P is a regular 2x2 matrix and M'{p) = Mi{p)+M2{p) is the sum of a diagonal matrix A^i and a Jordan- block 
matrix. M2{p) is chosen such that [Mi{p),M2{p)] = 0. M2{p) is nilpotent {{M2{p))'^ = 0). Thus, 
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which entries are: 
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The matrix Pip) is regular idetPip) ^ 0) if Mip) is regular. This decomposition leads to the form factors (p ^ 0) 
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Notice that X+{p) — A-(p) = \{p) = ^ imply the foUowing relations on the reaction-rates: 

E (■ 

9'=0,1,2 

E (. 



p/32p/3'l 
lO-"- 20 



p/31p/3'2\ ^ 
^ 00 ^ 00 ' ^ 



p/32p/3'l 
^ 10^ 00 



p/32p/3'l 
^ 00 ^ 20 



2/3pl;3' 

02 



pl/3p2/3'\ ^ 
^ 00 ^ 00 ' 



/3,/3'=0,l,2 



■ p2/3pl/3' 
^ 00^ 02 



pl/3p2/3' 
^ 00^ 01 



/3,;3'=0,1,2 



pippp z I ppipp z I p^lp2;3' . /pip . pPi\/p 

^20^01 "'"-'^02^01 "'"-'^02^10 ~r l_i 00 T~ ^ 00 



-^/3'2- 



< 



E/pl/3p2/3' 
20^ 10 

0,l3'=O,l,2 

00 + ^ 00 00 + ^ 00 J + 20 + ^ 02 Jl^ 00 +^00 J 

/3,/3'=0,l,2 

E/pl/3p/3'2 pl/3p/3'2 p2/3p/3'l p/31p2/3'\ _ /pl/3p/3'2 p/32pl/3 p2/3p/3'l p2/3p/3'A 

\^^02^10 T"^00^00 T'^OO^OO ^^00^01 j~ 2-^ \^^00^10 ^ 00^ 02 ^ 01^ 20 "'"^00 "^00 j 

■'=0,1,2 /3,/3'=0,l,2 

Erpl/3p/3'l , p2/3p/3'2 _ ,p2/3pl/3' _ ^pl/3p2/3' ,p2/3pl/3' ,pl/3p2/3' _ /pl/3 _ p2/3wp/3'l _ p/3'2N 
I -'^01-'^ 10 T"^02^20 ^^01^02 ^^00^00 ^-'^ 00 ^ 02 "T^^OO^Ol 1^00 ^ 00 00 ^ 00 / 

/3,/3'=0,l,2 

Erp2/3p/3'l pl/3p/3'2 _ /pl/3 _ p2/3wpl/3' _ p2/3' n _ /p/31 _ p/32wp/3'l _ p/3'2N 
I ^02^^10 "'"^01^20 ^ooA^^oi ^02/1 00 ^ooA^^io ^ 20 J 

/3,/3'=0,l,2 



/3,/3'=0,l,2 



(38) 



These are necessary conditions (but not sufficient) for the matrix M.(j)) to be non-diagonalizable. This means in turn 
that it is sufficient (but not necessary) that one of relations (38) be violated for the matrix AA{p) to be diagonalizable. 
If M2,i{p) = 0, the matrix AA{p) is triangular, i.e., 



A\ + Cl = A\ = Cl = Q, 



which in terms of reaction-rates imply 



E (■ 



,2/3 
10 



/3=0,1,2 

^2/3 _ p/32 



p/32 _ p2/3 _ p/32\ /p2/3 

^01 ^00 ^ m j ^ 2^ \^ 01 

/3=0,1,2 
n2/3 _ p/32 _ p/32 _ p2/3 



p2/3, _ 
^ 00 I ~ 



E 

/3=0,1,2 



r/32 _ p, 

^10 ^00 



(32 



(39) 



(40) 



So for example, if Tqq ~ Tq^ = T^^ — T^^ ~ ~ Fj^q , the relations (40) are fulfilled. 

If the matrix M. is diagonal, i.e., M^2,iij>) — ■A4i.2(p) ~ 0, and in addition to (|39|-^0|), wc have 



E (■ 

/3=0,1,2 



a/3 

20 



-^/31 
02 



pl/3 _ p/31 
^00 ^00 



E 

/3=0,1,2 



.1/3 
02 



-.1/3 
00 



E (rfo 



p2;3 _ ppz 

^10 — ^ 01 



/32 



As an example, relations (|39(-[4lD are fulfilled if one has T, 

p/32 

^ 10 



1/3 
00 



./31 
00 



.1/3 
02 



/3=0,1,2 



./31 
20 



.1/3 
20 



./31 
00 



and r, 



2/3 
00 



./32 
00 



(41) 



.2/3 
01 



it follows from this discussion that when the reaction-rates, in addition to the solubility constraints ([T^), also violate 
conditions (|3^) and one of the relations ( p8| ) which are sufficient but not necessary, the first case applies. When, 
in addition to (|l9|), the relations (p8|), are fulfilled (recall that (|3^) are necessary but not sufficient constraints) and 
the conditions ( p9[ ) are violated, then the second case applies. When reaction-rates satisfy ( p9| - ^0|) in addition to the 
relation (19), then the third case (see below) applies. 

Similarly, when reaction-rates satisfy (^9| - |4l| ) in addition to the relation (|l^, then the fourth case (see below) 
applies. 
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3) If A^2,i(p) = 0, the matrix Ai{p) is triangular, thus the eigenvalues of Ai{p) are A4i^i{p) and Ai2.2{p)- We have 
already discussed the physical implication of this case (see equation (|4^) above), and we have (p 7^ 0) 



which leads to {p^O) 



^ - Mi.iip) 



^M{p)t ^ gA^i,i(p)t ( 1 Mi,2{p)t 



S^{P.^) 



«>(0) 



uj - Mi,i{p) 



(42) 



(43) 



4) If both M.2i{p) = A^i 2(p) = 0, the matrix M{p) is already diagonal (see (EO)), and the form factors read 
(P^O) 



'^0 IP' ^) - uj-Mi,i(p) ' ' ~ ^-M2,2{P) 



(44) 



In the sequel we focus on the case where the matrix M.{p) is diagonal and provide explicit expressions in real 
space and time for the density. This is equivalent to imposing the six supplementary conditions characterizing Vi 
(see (HHl])), B'l + = = Dl = A\+ C| = Al = C\ = Let us first compute the density for the case where 
translation invariance is broken, i.e., the manifold Vi D V2 of dimension 72 — (16 + 6). With the above constraints, 
the densities obey the following equations of motion. 



-{nt) = (Ag + Co")d+ (n;l)(t)(A? + Q)d + 5] (A^(n;^+,.)(t) + C?(n;^_,„)(t)) 



(45) 



j^{nl) = d{Al + Cl) + {nl){t){B\ + D\)d + ^ (i?2^n^.+e= ) W + ^?("™-e=) W) (46) 

a 

In order to discuss the solutions of these equations, it is convenient to define the following quantities: ^.a = a/t4"j Ms = 

V 2 

Ca= J~^Cf et Cb = \/BlD\. Furthermore, we introduce Ba = 2(AJ + Cf ) <0, Bb = 2{B\ + D^) < 0, 

7A = + + + Q < 0, 7i3 = + + + < and = In^j, j = A,B. A site on the hypercube is 
denoted by m with d components: m = (m-i, . . . , m^, . . . , m^). 

Exploiting the well known properties of Bessel functions (j € (A, -B)), we arrive at 

(n^J(i)=p,(oo)(l-e-''l^^-|') + e-^l^l* ^ «,)(0) ^r^"-^— (47) 

7n — (?n^ ,. . . ,m^) i—l...d 

The steady densities are respectively, ^^^(oo) — 2\'yA\ ' ^^"^ Pb{oo) — 2\-tB\ ' ^^^^^^ investigate three different 
cases 

i) The (at t — 0) initial density (for each species) is given by pj , 

«')(* = 0) = Pj(0)<5™,„, (48) 
The asymptotic behavior is then {rrii 1 and uj — L"^ /4\Cj\t) 

«>(i)~P.(oo)+V^^ (49) 

where ipl,-^ are known functions. Further, we assume |7j| > 0, since for |7j| = 0, equation ( p^ ) tells us that (n^„)(t) = 
Pj (oo) = Pj{0) — cste,\ft, where 

e, =mm(d|7,|,d(^- |Q|(2 + eV2))) (50) 
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,ife, = d|7,|>o 

d/2 ,ife,=d(^-|Q|(2 + e2/2)) > 



For = and 7 7^ 0, the density decays algebraically, i.e., ~ Pj{oo) + '^^'^ . 

ii) Initially, the particles (for each species) are confined in some region of space, i.e., 

1^ , if < TOj < L/2 



«)(i = 0) = "0 ' " " - - ^1 " (52) 
^ ^ ^ ' , otherwise 



iiL/2 <m,<L 
otherwise 



We have then 



<(t)~p,(oo) + e-«^M^^;,+ 



where Oj has been defined in 



1) and for r = m — m', \r\ ^ 1, jCt < 00: 



, if = d|7j | > 

1 ,ife, =d(^-Q(2 + e2/2)) 



(53) 



(54) 



(55) 



If 9 = and 7 7^ 0, the density decays as a power law, i.e., {nl^){t) ^ Pj{oo) + {^m.i ^ — 
iii) The initial distribution of particles is assumed to be non- uniform, i.e.. 



«)(i = 0) = ^ "° 



if TO = 
if Iml > 



With (56), we arrive at (m^ > 1 and-Uj = L^/4:Cjt): 
Cj{2+ y)), decays as (p9|), with 



When = ^(1^^ 



E^f ,ifO<a, <1 
i , if a, > 1 



(56) 



(57) 



(58) 



Again, as Qj — 0, equation (p8|) holds. 



Notice the cross-over at 



= a 1, where the density decays as 



Pj(oo) + 



^_rf(^_|C,|(2+J-))t 



ln(4Mj|Cj|f) 
2-^4ir|C3|t 



By contrast, when Qj = rf|7j| > 0, the density behaves as {nl^){t) 



Pj{oo) + ■i/'^e"''!'''^!*. Here we have restricted our attention to the case where < a; < 1, Vi,0 < ai < 1, while in 
general we could consider different regimes in the different directions (for example, < ai < 1, ai = 1 and ai > 1 ). 
The corresponding asymptotic behavior follows as above. 

To conclude this section let us focus on the manifold V2, where the density is translationally invariant. On account 
of (^), we obtain two coupled linear differential equations which are easily integrated. The result is 



pj{t) = p,(oo) + (p,(0) - p,(oo))e-'^l^^- 



(59) 



The above solution allows us to solve for the correlation functions on V2{d). These in turn will be useful to solve 
perturbatively the problem on the manifold Vtransi-invid). 

Let us now discuss the relationship between our results and the solution of some models solved exactly in d > 1 
H earher. 
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In reference jO] , Clement et al., solved exactly (the fast adsorption rates version) of the Fichthorn, Gulari, Ziff 
(FGZ) model [112| introduced to describe the conversion of CO and O to CO2 on platinum substrates. The model 
solved in describes the dynamics of classical stochastic particles of two-species (called A and B) with hard-core 
constraint, and so our approach applies to this model. However, as the system solved in is a two-states model 
(in the stochastic variable is a "spin variable" Zj — ±1, +1 corresponding to an A particle and — 1 to a i? 
particle at site j. In this model there are no vacancies), it is not in the class of {three- states) models which we 
specifically study there. It is however possible to recover previous results |0 considering the system in the 
framework of the two-states analog of our approach |7|J|] . To do so we relabel Clement et al. B particles by vacancies 
symbolized by (the A particles are symbolized by 1), thus reactions occurring in are described by the rates: 

r}? = r?i = Til = ni = p/4; ni = rj? = i/4d; = t\1 = tI\ = r"? = p/4 + i/4d; r}} = r™ = -p/2 and 

FJq = Fq} = + 3/4d). For such systems (where 5=1), the 2s^ = 2 solubility constraints, which are the analogs 

of dll), read F™ + + F^J + FjJ = F?? + F?} + Fj? + Fj} and F™ + FgJ + Tj? + = T?? + Fj? + T% + F}J. 

These constraints are fulfilled for the previous choice of rates (similar to the choice of |Ol) and thus the dynamics 
of the system is soluble in arbitrary dimensions, i.e. the equations of motion of the correlation functions are closed 
and obtained as in ( ^JitIJi^ ). As an example, for the density we have the following equation of motion Q: ji{nj) — 
Bi J2±a ~ ^j)) ~ P{''T-j)j whcrc Uj = Uj — ^ and Bi = F^J + FjJ — Fqq — FjJ — l/2d are the same quantities 

defined in Noting that in language of ||ll| = (zj) ~ 2{nj), we recover the result of Q: Jj7j = -^^Ij —Plj, 
where /S^j = X]±Q(7i+ct — Ij)- Similarly we can reproduce the (closed) equations of motion for higher correlation 
functions. Saturation phenomena as in ||l^ should also occur in the class of three states models. However the 
analytical treatment would be more complex than in | [Tl[ | two-states models. 

Let us sketch the strategy which one should follow to treat saturation in the models considered here. For trans- 
lationally invariant systems one should solve a (linear) differential-difference systems of coupled equations describing 
equations of motion of correlation functions {n\^n\){t), j £ (A, B) paying due attention to the boundary terms m ~ I 
and |m — Z| = 1 (see also section V). This system is solved in Fourier space and involves a general 3x3 matrix with 
non-constant entries. One should, as it has been done for the density, carefully discuss the properties of this matrix, 
which is a technical matter. In fact such a study should be carried out for some specific model. 

A further two-states model which can be solved exactly in d > 1 is the Voter model (see e.g. [|l^) described 
by the reactions-rates: F = Fq^ = F^g = Fj} = Fj} > 0. Since this model fulfill the previous "two-states" 
solubility constraints, it is soluble in arbitrary dimensions. With the two-states analog of (^ 17p^) we obtain 
the (closed) equations of motion of the correlation functions. As an example, for the density, we have [^ : 
^(n,) - 2F^„ ((n,+,) + (n,_,) - 2(n,)). 

Another important model which has been studied rigorously in dimensions d > 1 is the (irreversible) reaction 
A-\- B ^ % -\-%. For this model, Bramson and Lebowitz |Q obtained, rigorously, upper and lower bounds for the long 
time behavior of densities. However, systems considered in these works, 10, allow the multiple occupancy of a site 
by particles of the same species. Later, Belitsky |^ generalized the study of jlj] to the case of hard-core particles 
reacting according to A + B^% + ^ (with rates I^§ ^'^21 = r(= 1)) and A + ^ + A; B + % ^ % + B (with 
rates T\l ~ Fq" = Fjo = Fq2 = 1). He obtained rigorously an upper bound for long-time behavior of the density 
[for an uncorrelated initial state with equal species densities: /9a(0) = Pb{^) < 1/2]: Ve > 0, 3T(e) < 00, so that 
for t > T(e), pt^oo(i) < for d < 4 and pt^oc{t) < C*t~^, for d > 4; where C* is a positive constant. We 

can now wonder whether such a model can be dealt with in our approach. As the model considered by Belitsky is 
a three-states model, the equations of motion of correlation functions are given by (|l4|,[17[[lq) with A2 = Cf = B\ 



-"2 

D\ = ~A1 = -B\ = -Cf = -Z?^ ^ I and Al = BI = B^ = Dl ^ ^ Al ^ ^ ^ C\ ^ = Q. 

Unfortunately the solubility constraints (|l9|) are not fulfilled for such a model (consider e.g. the fourth constraint: 
A2 Bf -\- Aq = F2]' -I- F21 -I- Fgf ^ F = 0, but in this model F = 1 ) and the equations of motion give rise to an open 
hierarchy which cannot be solved. 

Results in d > 1 have also been obtained by approximate methods (mean-field theories), and/or by scaling and 
heuristic arguments (see e.g. ]l| and references therein). As an illustration of these studies let us consider the work 
of Toussaint and Wilczek. In pS, a system of two species A and B reacting according toA + i3^0-|-0, in addition 
to their diffusive motion is studied numerically and an approximate method for calculating the densities at long- 
time is proposed (for system with equal densities : pa(0) = p_b(0)). Approximate results Q predict p{t) ^ t^'^^^ 



rigorous results in the one-dimensional case, but in 
|13| is a continuum macroscopic approach and cannot 



in agreement with Bramson, Lebowitz ||1J] and Belitsky 
disagreement in higher dimensions d > 1. The approach in 
take into account the hard-core constraint of the particles . In addition, it takes into account of fiuctuations in an 
approximate and uncontrolled way. It is therefore difficult to compare their method with the microscopic exact results 
presented here. 

In sum we have seen that the two-states formulation [Q of the method dicussed here allows to recover some previous 
exact results in arbitrary dimensions pT| , p^ for the stochastic models of hard-core particles. Our approach applied 
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on three-states models (of hard-core particles) is in a sense complementary to the rigorous results of jlj] and is useful 
to describe exactly physical models such as a three-states growth model [|lO| . 



IV. NON-INSTANTANEOUS CORRELATION FUNCTIONS 



To our knowledge, the only exact computations of two-time correlation functions in non-equilibrium statistical 
mechanics that are available are those for single species models, in particular for one-dimensional models which can 
be mapped onto free fermion [ p8| , p"9p^ ] and other related |2^,|l^ models . These exact results are useful, as starting 
points for perturbative calculations or for checking numerical computations. This section is devoted to the study of 
the 50-parametric manifold Vi . We are interested in the density-density correlation functions 

{nl^it)nim EE {x\nl^e-"'ni\Pm (xKe-^*|F'(0)), ^,J e (A, B) (60) 

where \P{0)) denotes the initial state of the system. From the above we see that the evaluation for the correlation 
function with respect to the initial state |^'(0)) is equivalent to computation of the density of particles of species i at 
site m for a sytem in an initial state described by |P'(0)) = nj\P{0)). 

We now distinguish the case where correlations are absent in the initial (with broken translation invariance) state 
from that where they are present. 



A. Non-instantaneous two-point correlation functions for uncorrelated initial states 

In this subsection we assume uncorrelated initial states with a random and translationally invariant distribution of 
particles of type A (density Pa(0)) and of type B (density pb{0))- Therefore in our notations \P{0)) becomes 



So that we have 



I - PAiO) - pb{0) 

\pm = I pa{o) 



{ni{Q)nf{Q)) = ^^(0)5™,, + pa(0)'(1 - 5^,i) 
{n^iO)nfiO)) = pB{0)6rn.i + Pb{0)'{1 - <5„u) 
{ni{0)nfm = {ng{0)nfm = pa(0)pb(0)(1 - S.^j) 



(61) 



(62) 



We begin this section by computing explicitely the Fourier-Laplace transform of the two-points non-instantaneous 
correlation hmctions, i.e., the dynamic form factors measured in the light scattering experiments |r^Jl9t] 



1 

^ / ^„ J I' 



m'— m— / 



m' —rn — l 



dte-^P-^'-'{nln'=m-iitKmAhj) e {A,B) 



(63) 



Using the results of the previous section and assuming regularity and diagonalizability of A4{p), we find for the 
dynamic form factors 



1 



X.{p)-Mis{p) \+{p)-MiAp) 



X-{p) - X+{p) [ uj-X+{p) uj-X-{p) 

•Mi,2(p) (PA(0)pa(0))(Vo~l) 
[lo - X+{p)){io - X^ip)) 

[(7- - MuimiA^o + CS)d - Mi,2(0)(Ag + C'o)d] 

7- - 7+ 

[(7- - MuimiA^o + CS)d - Mi,2m< + C!^)d] 

7- ~ 7+ 



{{p\{0) + PAmSp,o - pa{0)) 



1 



1 

uj{uj — 7-) 



(64) 
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and 



A_(p)-A+(p) 



-X-ip) 



PA{Q)pBmSp,o ~ 1) 



+ 



Mi,2{p) ((p|(0) + PsmSpfi - PBjO)) 
{uj - X+{p)){co - X.{p)) 



7- -7+ " . • ■ ■ " . cj(cj-7+) 

[(7- - Xi,i(0))(AS + CS)d^M^,2m< + C'„)d] --^ ^ 

7- -7+ uj(uj-j-) 



X^{p) - Mi.iip) X+{p)-MiAp) 



u- X- {p) 



^ - A+(p) 



{X+{p) - MiAvm-ip) - Misip)) 



Mia{p){Lo - X+{p)){cu - X-{p)) 



1- - 7+ 



7- - 74 



(7- --Mi.iW) 

Xl.2(0) 

(7+->(i4(0)) 

Xl,2(0) 



{A^^ + C^,)d-{Al + Cl)d 



p^(0)pB(0)(<5p,o - 1) 

7+ -Xi.i(O) 



7_ -Xi,i(0) 
a;(a> — 7_) 



(65) 



(66) 



5f^(p,c^)-p^(0)pB(0)(Vo-l) 



A_(p)~A^i,i(p) A+(p)-Ali,i(p) 



w - A„ (p) 



(A+(p)-.Mi,i(p))(A-(p)-A^i,i(p)) 
-Mi,2(p)(c^-A+(p))(c.-A_(p)) 



((pi(0) + pA(0))<5p^o-PA(0)) 



7- - 7+ 

1- - 7+ 



^%y^ A- + Cg)d-(A^ + Cg)d 



7+-.Mia(0) 
w(w - 7+) 

7_ -7^1,1(0) 
cj(aj — 7-) 



(67) 



Again the poles of the dynamic form factors give the relaxational eigenvalues. As in the previous section, we could 
also compute the correlation functions in the case where M.{p) is non-diagonalizable, triangular or already diagonal, 
but for brevity's sake we prefer to focus here on the non-instantaneous correlation functions on the 50— parameterss 
manifold Vi. 

With help of equations ( p5[]47| ), we obtain the non-instantaneous two-point correlation functions on the manifold 
Vi as 

("™(iK(0))=PA(oo)(l-e-'^l'^-l*)+pi(0)e-''l^-l* + (p^(0)-p^(0)) J] A^:,"°~'"e-^/„,_i J2a^i) (68) 

a = l...d 



'Mnfm = Pa(oo)(1 - e-''!^-!*) -t- pa{Q)pb{Q) 



^-d\^A\t 



n 

a=l...d 



e 2 / 



ra^j —I, 



X2CAt) 



(69) 



(n^(i)nf (0)) = ps(oo)(l - e-'^'^^l*) -f p|(0)e-'^l^«l' + (ps(0) - pUO)) JJ <-'°e'^/„.„-, j2Cst) (70) 



(n^ritWm = Pb(oo)(1 - e-'^l^-l*) + pa{0)pb{Q) 



-<i|7Bl* . 



^i...d 



(71) 



with the notations, r = X^Q^aCa, el = eiX^Q^Q and Tq, = mc — la — crL. We are interested in the asymptotic 
behavior {\Cj\t 1, j G {A, B) and Wj — L^/4\Cj\t < oo) of the above correlation functions in two regimes: 
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1 /2 

i) when \m - l\ = \r\ = {Y,a=i...d '^D ~ i > 1, in this case cr = r/L = 

ii) when \m — l\ = \r\ <^ 1, in this case a — r/L — 0{1/L). 

It is worth noting that the autocorrelation functions (where |m — Z| — 0) are obtained in the second regimes (ii). 
With the above, we finally arrive at 



(PA(0)~-Pi(0))e-'^"'"" 
{A7r\CA\tr/^ 



(72) 



{ni{t)nf{0)) = pa{oo){1 - g-'^l'^^l*) + pA(0)ps(0)e-''f^^l* + exp ( 2da^7.A 



X PA{0)pBiO) 



1 



(47r|CA|i)''/2 



C'(l/i'') 



|C^|(2+4/2))t 



(73) 



"(PB(0)-p|(0))e-''-V 



exp 2da ub 



2|CB|i 



|CB|(2+4/2))t 



(47r|C^|t)'^/2 



0(1/*'^) 



(74) 



(n^,(t)ni4(0)) = pb(oo)(1 - e-"^!^-!*) + p^(0)pB(0)e-'^l^-l* + exp ( 2da2uB - (LJe}£^\ ,-d(^-|c.|(2+4/2))t 



2|Ci3|i 



X pa(0)pb(0) 



(47r|CB|t)''/2 



©(l/i'^) 



(75) 



In the regime \Bi\ = 1^^1(2 + ef/2d), i E {A,B), the two-point correlation functions decay as 



exp (^2da^Ui 



i?-e,\c,\ty 

2|Ci|t 



(z,j)e (A,s) 



(76) 



Note the non-trivial dependence on dimensionality and the drift term in exp (^2da^Ui — ^^~2\c^\t*'^ 
remark that this is consistent with the result obtained in one dimension for free fermions ||l9|. If e^ = 0, than there is 
no drift: exp (2da^Ui — ^^—^f^r^^j = 1- 



, {e^ 7^ 0). We 



B. Non-instantaneous two-point correlation functions on Vi: correlated initial states 



Let us consider correlated initial states described by a distribution having the following properties: 

i) when dist{m — I) > Q: 

«(0)n^(0)) = JC,j W {I- <^,.„,o)|ra|-^--, Ay > > ,r, = |to, - € {A,B), (77) 

a=l,...,d 

ii) when m = I, 

(KMnim = {KM)s^.J - (W(o))*,, = pms..j (78) 

The initial distribution in this subsection has been chosen in a special form, namely in such a way that computations 
can be carried out explicitly to the end. There is a priori no physical justification for such a choice, which has already 
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been considered in pl[ | for a single-species reaction-diffusion system. However and most importantly, our goal here is 
to investigate the dependence of the asymptotics on the initial correlations. 

We remark that in one dimension the initial state (^;^8|) is translationally invariant, 

{nl^{0)nj{0)) = (n|,|^|„_;|(0)n^o(0)) = /C,,(l - S^ri,o)\r\-^^^ + p^{0)6.,,6^rl,o (79) 



while in higher dimensions (d > 2, see 77) translational invariance is broken. This state of affairs has lead us to 
distinguish in the discussions the one one dimensional case from its higher dimensional counterparts. 

i) We begin with one dimension {d = 1), where r — = m — l. Because of the translational invariance in the initial 
state, we expect the non- instantaneous correlation function to depend on r ~ m — I, indeed, 

{ni{t)nf{0)) = {n^{t)n^m = Pa(oo)(1 - £-1^-1*) + p^(0)e-^/i!4/r(2C^i) 

+ ICaa J2 t^7'''yr^-'^e-^Ir-r'{2CAt) (80) 

{nt{t)nfm = p^(oo)(l - e-l^-l*) +/Cab ^ Ma^'^V' I "^""e"^ (2^ (81) 

{nl{t)nfm - {n^{t)n^m - Ps(oo)(l - £-1^-1*) + ps(0)e-^Ai's/,.(2Cst) 

+ ICbb /i^-''^>r^""e-^/,_,.(2CBi) (82) 

{nf^{i)ntm - Pb{^){1 - e-l^-l*) +/Cba ^ A^B"''Vr^""e-^/._..(2CBt) (83) 

Notice that when /i^.s 7^ 0, then (n^'-^(t)n^'^(O)) ^ (n^^^(i)n^''^(0)) because of the drift which is due to an 
asymmetric Markov generator H . Such a behaviour has been observed in single-species one-dimensional free-fcrmionic 
models fllJll. 

ii) In higher dimensions {d > 2), the initial state is no longer translationally invariant and with ( [77| , [78| ) and (|4^), 
we find 



d\BA\t 

a=l...d 



+ ICaa n Mr"'"1<"^or^^-^e-^/„„_„^(2CAi) (84) 

(m[^h,...,m'^^ld) a=l...d 

+ ICab Y n A'"°""°K-^ar^^-e-^/™„-™'j2C^t) (85) 

{■m{^li,...,m'^^ld) a=l...d 



(«™(0"f(0)) =PB(oo)(l-e-'^l^-l*)-|-ps(0)e-^ Y[ (2CBt) 

i...d 



+ JCbb J2 n A^s°""°K-^ar'^^"e-^/™„_,„.j2CBi) (86) 

(m[^h,...,m'^^ld) a=l...d 

{n^{t)nfiO))^PBi^)il-e-'\''-\') 

+ K.BA Y n A*B°"''"1<-^ar''^"e-^/„„_„, (2Csi) (87) 

(m[^h,...,m'^^la) a=l...d 
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We see that in higher dimensions, because of the broken symmetry of the initial state, the non-instantaneous correlation 
functions no longer depends on = TOq — la- We can study the asymptotic behavior of these non-instantaneous 
correlation functions in an unified way (including both d = 1 and d > 2), namely for — — la, with \ra\ = 
{iTT-a — la\ ~ 3> 1, with — OaL and Ui — /A\Ci\t < oo, in the regime where |Ci|t, r ^ 1, G {A,B), the 
correlation functions are given by 



2\CM 



n 



= l...d 



©(i-^'') , < A" < 1 (88) 



where we use r = J2a "^a^a and el = Moreover, 



«(iK(0))=p.(oo)(l- 



^-'i|7^|t^ 



When = 1, a cross-over takes place and logarithmic corrections arise, namely. 



(89) 



A" = 1 



(90) 



Therefore, when the spatial correlations are important, i.e., A^- < 1, correlation hmctions (at — \Ci\{2 + ef /2d)) 
0) decay as 



("LW"l(o)) 



(91) 



(47r|a|t)S.^S/2 

On the contrary, weak spatial initial correlations do not affect the long time behavior of correlation functions, since 

•{^Y.a-lu.~{'-^) 



exp I 



A« > 1 



(92) 



The marginal case A"- = 1 has logarithmic corrections 



«WW'(o)) 



(2 Eg <^ln^ - {^^0^?) (ln4|Q|t)° 
(47r|a|t)^/2 



A" = 1 



(93) 



Notice the drift which occurs for e 7^ and the effect of dimensionality. The fact that initially the state is 
translationally-invariant {d = 1) gives rise to the same asymptotic behavior {t~'^/'^) as for the non-translationally- 
invariant system in higher dimensions {d> 2). 

V. INSTANTANEOUS TWO-POINT CORRELATION FUNCTIONS ON THE MANIFOLD Vo 



We now pass to the computation of the two-point correlation function on the translation invariant manifold V2 (d) 
(pT|). From (^ 18), the evolution equations of correlation functions follow. We shall discuss both cases, when the 
initial state is correlated and when it is uncorrelated. We shall evaluate 
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qAA 










qBB 
=^\r\ — \n- 




{n^ng){t) = 






=^\r\ — \n- 


-rn\it)= 






(94) 



with the boundary conditions at r = (for the densities see (^ 

G^^o(t) ^ PAit) ; C/ffo(0 EE pB{t) ; G^^it) ^ 



(95) 



Notice that the point \r\ — 1 must be dealt with care. Further, it wiU be convenient to distinguish the one dimensional 
problem from that in d > 2. In this section, in addition to t he defi nition of th e Appendix A, we also introduce some 
specific notations and abbreviations (see Appendix B, (143), ( 144 ) and ( 145 )). Let us start the discussion with one 
spatial dimension. 



A. One-dimensional instantaneous two-point correlation functions on V2{d = 1) 



On account of the above remarks and (|17|), the equations of motion for the correlation functions read (|r| > 1): 

J^Q^^{t) = BAQ^^{t)+CA {g^Aii)+QrA{t))+AAPA{t) (96) 

while for \r\ — 1, we have (|l|), 

J^Gt^{i) = (Gi + BA/2)gf^{t) + CAG2it) + + (Ff + + ^A/2)pA(i) + {F^ + Ft)pBit) (97) 



For r = 0, we recover 



rAA(.-. d Aa , fBA + 2CA 



PA{t) 



(98) 



The solution of the above set of coupled differential equations (|9q-p8D can be expressed in terms of the modified Bessel 
hmctions Iv{z) (See appendix C), and for \r\ = |n — to| > we have: 

g^'^it) ~ {pA{t)f - -(pA(0)e-l^^l*)2 + {pa{0) + ^^£±^1±^ ] e-l^^l*/,(2CA0 + Qr'^i^y-^""''^' Ir-r-{2CAt) 



+ + 



A . T^2,o\Ba\ 



Ca 



Ca 



H,{\Ba\~\ia\) 



Ca 



e-lT-^l* / dTe-(l^l-lT-*l)"/,(2C^T) 







+ 



' vU\Ba\~\ib\) 
Ca 



j e-l''^!* dTe-(l^^l^l''«l)^/^(2CAT) 



similarly, we find 



{G't-BAl2) d 



Ca dt 



-Ir{2CA{t - t')) + 2CAlr{2CA{t - t')) 



(99) 



Q?^{t) - {PB{t)f = -(ps(0)e-l^-l*)2 + pb{0) + 



Cb 



\*Ir{2CBt) + J2 Q?''iO)e-\''^^'lr-r'i2CBt) 



Cr 



^ dTe-l^«l^/,(2CBr)+ {vl^ 



T^U\Bb\-\ib\) 



j e"l''«l*^ dTe-(l^«|-|^«l)^/r(2Ci3T) 



VU\Bb\~\ia\) 
Cb 



rdt'e-l^-l(*-*')5f^(i') 

-'0 



-|7A|t 



^^g-(|BB|-|7^l)rj^(2Ci3r) 



{GI-Bb/2) d 



Cb dt 



-Ir{2CB{t - t')) + 2CBlr{2CB{t - t')) 



(100) 
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and we also obtain 



/t)AB I jjAB I jjAB \ 



Cab I 

-D^filBABl-bAl) 



2Cab f dt'e-\^'^^\^'-''^gt^{t')Ir{2CAB{t^t')) (101) 







To study the asymptotic behavior of these expressions, we shaU distinguish two regimes: 

i) long time , i.e., \Cj\t S> 1 and large distances , i.e., r ^ 1 with the ratios j^-j^ < oo and = L? /A\Ci\t < oo 
hold finite. 

ii) long time , i.e., \Cj\t ^ 1 and finite distances, r ^ i ^ oo with j^-^ ^ 1 and ui = L^/4|C;|i < oo. 
In order to investigate the effect of initial correlations on the dynamics, we consider 

^M>o(0) = P»(0)P,(0)(1 - \Ki\{sign{Ci)Yr-^'), vi > 0, (102) 

where {I e {AA, BB, AB), i,jeA,B). 

Such a choice has been made for the one-dimensional single-species symmetric A + A ^ ^ + ^ process . In this 
subsection, we want to proceed with a systematic study of instantaneous correlations, on the manifold V2{d = 1), 



with the choice (102) for the initial state. 

Below we shall use the incomplete gamma function, T^v, u) = dxe^'^x'^^^ , as well as the Riemann zeta functions 
C,{y) = 'Ylik>i^~'^ ^ ^- ^'^'^ notational simplicity, we write v instead of vi. The results for the asymptotics are 
summarized" as follows 

1. For \Bi \ > 2\Ci\, the decay of correlations is exponential. With the definition 



^,,=rmn{\j,\,\B,\~2\q\) (103) 



we have 



l.a.Foripjj^\^A\^\Bj\-2\Cj 



Qi-Ht) - g^^^oo) e-\''^^' (104) 



l.b For ipjj = \Bj\ ~ 2\Cj\ > 0, and g\r\>(){0) + 

gi'{t) - g«'(oo) ~ p,(0)2e-(l^^l-2|c.l)t (105) 
Note that in the case where t/r(0) — 0,Vr, we have 

g^(^)~g^{^-\ Justus,-.,,, . (106) 

^ ^- ' SI 7- 1 



18 



where Q\{uj, a) is a function explicitely determined by the processes which occur in the system under consideration. 

l.c For \Bj\ ^ 2|Cj|, Yij'=/.i\ = \Bj\- 2|Cj|, |7j| > 0, and V\.^ + 2{signCj)V^^^^ ^ 0, we have 

f e-lT^I* if = |7,| 

«fW-«'(oo)~|,Hv„,.;^it,, = i7.«i 

l.d For \Bj\ ^ 2|Cj|, |7j| = \Bj\ - 2|Cj|, |7j/^j-| > 0, and V{^^ + 2(si5nCj)I?^ ^ 0, we have 



e^' W - Q'A^) ~ (108) 

The correlation functions Q^^ (t) have to be discussed separately. With the definition 

^AB = mm(|7A|, |7i3|, \Bab\ - 2\Cab\) (109) 
l.e For (fiAB = |7a| \Bj\ - 2\Cj\, we have 

gf''{t)-gf''{oo)r^e-^''^\' (110) 

l.f For (fiAB = |7b| 7^ - 2|Cab|, we have 

gf''{t)-g^''{oo)^e-\-"'^* (111) 

l.g For .fiAB = \Bab\ - 2\Cab\ > and g^^g ^ 

a,^^(0 - e?^^^(oo) ~ pA(0)ps(0)e-(l^^«l-2|C^B|)t (112) 
Notice again that if g^'^iO) = 0, Vr, we arrive at 

' Qi(u.i,a)e-'^''ABl-^lCABl)t 



grit)-gr{oo)^< 



, if r < L 

Q,iuAB,T)e-''^"^B-'.\BABl-^\CABl)t ^ ^ ^ ^^^""^ 

^/4■K\CAB\t 



Where Qi{uj,a) is the same quantity as above. 

l.h For \Bab\ ^ 2|Cab|, |7a| = \Bab\ - 2\Cab\, and V^f + 2{signCAB)'Dif ^ 0, we have 

l.i For \Bab\ + 2|Cab|, |7s| = \Bab\ - 2|Cas|, and + 2{signC ab)'^^^ + 0, we have 



c^^m c^^roo^ fe-i-i*^ if^AB = l7A| 

Qr it)-Gr (-)-|e-l-l*^if^.B = |7B| 



(114) 



(115) 



2. For = 2\Ci\, the correlation functions decay algebraically. It is appropriate to distinguish, r <^ L and r ^ L. 
Again, a = r/L, I e {AA, BB, AB). 
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In the regime r <^ L, 



(4|Ci|t)-/2 

[2CM + (4<T"n,|C,|t)''-"'^"] 

(47r|C,|t)i/2 
In [4|Ci|«l(l-g-)t] 
(47r|C,|t)l/' 



if < 1/ < 1 

if i/ > 1 
if z/= 1, 



where the following auxiliary function has been used 

Ti (u, a, v) = ^ 



47r 



while for r 3> 1,7" = crL ^ L, we find 



with the auxiliary function 



(4|C,|t)-/2 

[(1+e— '"i)C('') + ((l-'^)/<^)(4o-^«i|C,|t)<'""'/'] 

(47r|C,|t)i/2 
ln(4|Ci|Mit) -a-^Mi 



if i/ = 1, 



J^2{UUCF, V) 



Notice that when the initial correlations are absent (i' = 0), 

' Qi(u,a)+0(l) 



^/iTT\Cl\t ' 



S\Ci\t 



if < 

if r < L and C; > 



1 and C, > 



(116) 



(117) 



(118) 



(119) 



(120) 



where Qi{u,a) has been defined previously. 

From the above we infer that th e ini t ial condit ions affect the asymptotic behavior of correlation functions only 
when the latter decay algebraically (HE, IIS, 120). Provided the initial correlations arc dominant, {0 < vi < 1), the 
critical exponent is renormalized, while for weak initial correlations, {i>i > 1), the exponent is independent of initial 
correlations, i.e., 1/2. The intermediate regime, vi = 1, exhibits logarithmic dependence consistent with a marginal 
behavior. 



B. Two-point instantaneous correlation functions on 1^2 (d) in arbitrary dimension 

This section is devoted to the computation of correlation functions in arbitrary dimensions (d > 2) on the manifold 

V2id). 

With the notations r = (ri, . . . , Tq,, . . . , r^), |r| — \/J2a (sometimes denoted by r, for notational simplicity) and 
= \J [toi ± 1)^ + 'Ylia^a' ''a'' solviug the multidimensional equations of motion of the correlation functions (see 
appendix D), we arrive at the following explicit forms 



= l...d 



r'=iO a=l...d \ 

' PU\Ba\-\ia\)^ 



VUBA\d\ 



Ca 



U'^ n 



-I^-^I*/^„(2Ct) 



+ d- 



Ca 



g-<i|7A|t / ^^^-d{\BA\-\lA\)r -Q 1^^{2Cat) 
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+ dV^^e-''\'rB\t( \BA\-\lB\ \ j\^^_a(\BM-\iB\)r TT 7^J2C^t) 

-|*diU(i')e-'^l^^l(*-*')(2C^d n ^-(2C4i-0) + ^^f^|7 n IrA2CA{t-t'))Y (121) 

^ \ c«=l...d a=l...d / 

we also have, 

fCBPA{0)+Vf.+Vf,+VJ' 



r'T^O c<=l...d V B J Jo a=l...d 

- j\t'g,{t')e-'\^'^\^'-''^ i'iCBd n ^-(2CB(i-0) + ^^f^|7 n /.„(2CB(t-0)), (122) 

^ \ a=l...d c«=l...d / 



and 



a^^l(,,,...,,,)l>o(*) - (PaWpbW) = -(pa(oo)ps(O) + p^(0)pb(oo) - pA(oo)pB(oo))e-l^-+^«l'^* 
-(p^(0)pB(0))e-<^(l--l+l--l)*+(^^^i±||J±^j n (e-l^-l'/.j2C^Bt)) 

+ E^m(0) n (e-l^-l*/..-.^.(2C^Bi)) + (2?5!^ + ^^^^) /dr [] (e-l^^-l%j2C^Br) 

a=l...d \ / ^ a=l...d 

V ' Cab ) Jo J,\^ 

- f dt'g,{t')e-'\^--\('-'')UcABd n J.. (2C^B(t - ^0) + ^ ~ It n IrA'iCAB{t-t'))\ (123) 

V c«=l...d a=l...d / 

Notice that these expression are vaHd, on Viiii), in arbitrary dimension and setting d = 1 we recover the one- 
dimensional expressions of the previous section. 

We assume here that the initial state is characterized by a random, translationally invariant, but uncorrelated initial 
distribution, i.e., fj[^|^(,(0) = pi(0)pj(0). 

The asymptotic behavior is obtained similarly to the one dimensional case when \Ci\t ^ 1 with ui = L"^ /A\Ci\t 
1. For \Bi\ > 2|C;|, the decay of the correlation function is exponential, 

<Pjj=min{\jj\,\Bj\-2\Cj\) (124) 

l.a If fjj = \ja\ 7^ \Bj\ — 2\Cj\, we have 

e^^(0-er'(oo)~e-'^'''''* (125) 
l.b If ipjj = \Bj\ - 2\Cj\ > |7j| and e^^o(O) 0, we have 
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QiHt) - QiKoo) ~ e-''(l^^l-2|'^^D* (126) 
Finally, note that provided ^^^0(0) = 0, we find 

I (47r\Cj\ty/^ It r » i 

where Q2{uj,aa) is a known function determined by the processes occuring in the model under consideration, 
l.c If \Bj\ ^ 2\Cj\, |7jv^j| = \Bj\ - 2\Cj\, |7j| > 0, and 2 + 2{signCj)V{ .2 ^ 0, we have 



{t) - gi\Qo) < 



e 



e-^l-r.v.l* int , if (^^j- = and d = 2 

, e-'*lT/^il*(47r|Cj-|i)i-<^/2 , if Lp^j = and d > 3 



l.d If \Bj\ ^ 2|Cj|, |7j| = - 2|Cj|, > 0, and Dj ^ + 2{signCj)V{ -^ 7^ 0, we have 



e-2|^^l*lnt ,ifd = 2 (129) 

_ e-<i|7.l«(47r|C,-|t)i-'^/2 , if rf> 3 

The functions g^'^H) requires a separate discussion. With the definition of ^pab 

ipAB = minQjAl, \7b\, \Bab\ - 2|Cas|) 
l.e If (fiAB = |7a| 7^ \Bj\ - 2\Cj\, we have 

g^'^it) - g^^'ioo) ~ e-'^lT^I* (130) 
l.f If ifAB = |7b| \Bab\ - 2|Cab|, we have 

Or^'it) - g^^'ioo) ~ e-'^l'^^l* (131) 
l.g If ipAB = \Bab\ - 2\Cab\ ^ |7a,b| > and ar#o(0) ^ , we have 

g^^it) - g^'^ioo) ~ e-'^(l^^^l-2|CAB|)t (132) 
Again, when ^^1^(0) = 0,Vr, we arrive at 



rABf,\ r^Br-. I (47r|CAB li / <^ ^ 

yr [t) - g^ (00) ~ < Q,(„AB,a„)e-''"AB-.i(IBABl-2|CABl)« 

t (47r|CAB nr^i 

Where Q2{uab, era) has been defined previously. 

l.h If ^ 2\Cab\, \ia\ = \Bab\ - 2\CabI and V^^f + 2{signCAB)Vlf ^ 0, we have 

" e-'^lT^I* ii^AB = \lB\ 

= -|7A|t / t 



(133) 



g^^{t) - g^ioo) ^ ' V '^^^ = ^ ^ ^ (134) 

g-2|7A|tini if (^^B = |7a| andd = 2 

[ e-'^M*{4Tr\CAB\t)^-'^/^ if ^as = |7a| and > 3 

l.i If \Bab\ ^ 2|Cas|, |7s| = I-BabI - 2\Cab\, and Df^f + 2(signCAB)^?^f ^ 0, we have 
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AB, 



-d\^A\t 

\lB\t 



e 

e I'^^l'lnt 



\lB\t 



-d/2\ 



if ipAB = \^a\ 

if ipAB = IT-bI and d = 1 
if (^AB = |7b| and d = 2 
if (^AB = |7b| and d > 3 



(135) 



2. For = 2|C;|, the decay of correlation functions is algebraic. We distinguish the regime r L, from that 
where r ^ L {a ^ r/L, I G {AA, BB, AB). 
In the limit r <^ L, 



if C; < 



8\Ci\t 



(i-n\Ci\t)'i/^ 
; {ui ^ ) 
(47r|C,|t)''/2 

with the definition of the following auxiliary functions: 
In the regime where r ^ 1, Tq, = cTai ^ L, wc find 



47r 



P.(0)Pj(0) 

sl'Cilt 



(47r|C,|t)''/= 



if C; < 

ifCi > 0, 



with the auxiliary functions defined as: 
Notice that setting ^^,.(0) = 0,Vr, leads to 



gUt)-gl{^) 



(47r|Cjt|)<i/2 

Q2(Ml,o-Q)e ^t.' 
(47r|Cji|)''/2 



if r < L 
if r > 1 



(136) 



(137) 



(138) 



(139) 



(140) 



where the function Q2iu,a) is as above. 

We see that for the uncorrelated initial state under consideration, the dimensionality of the problem has a non- 
trivial effect on the dynamics. In fact, in the critical regime (138, 14C), when d> 2 the correlation functions decay 
as instead of as in lower dimensions. We further remark that also in the mas sive case , the dime nsionality 

of the model can have particular non-trivial effects on the asymptotic regime (see e.g. 128, 12£, 134, 135 ). Let us 
conclude by noting that all results obtained in this subsection are compatible and in agreement with the previous 
one-dimensional results. 



VI. CONCLUSIONS 



In this technical paper, we have classified the solutions of the two-species bimolecular diffusion-limited reaction 
models and have been able to obtain exact and explicit results, namely : 

- The Fourier-Laplace transform of the density and of the non-instanateous two-point correlation functions (dynamic 
form factors) on a 56— parameters space, in arbitrary dimensions. 

- Exact computation, in arbitrary dimensions, of the density on a 50— parameters manifold for various initial 
conditions. 

- Exact computation of the non-instantaneous two-point correlation functions on a 50— parameters manifold for 
uncorrelated homogeneous, but random, initial states as well as for initially correlated states, in arbitrary dimensions. 
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- Exact results for the instantaneous two-point correlation functions on a translationally- invariant 31— parameters 
space manifold in arbitrary dimensions. 



Exploring the various classes of solutions for the one- and two-point correlation functions, we have seen in real space 
and time that there are essentially two regimes, a massive one and an algebraic one for the density and two-point 
correlation functions. For non-instantaneous correlation function we have pointed out that a drift can occur due to 
an asymmetry of the reaction-rates characterizing the stochastic Hamiltonian. We have also shown that when initial 
correlations are strong enough, the critical exponents in the asymptotic regime are renormalized while for weak initial 
correlations, the long-time behavior is insensitive to and independent of the initial state. 

Our approach applied on three-states models (of hard-core particles) is in a sense complementary to previous 
rigorous results (see the end of section III) and allows to study exactly physical models such as a three-states 
growth model ]10|. From our analysis, it follows (alternatively, using the Rauth-Hurwitz conditions and simple 
algebra) that in arbitrary dimensions, two-species models belonging to the class of soluble models discussed here, do 
not exhibit phase-transition, nor pattern formation. This has lead us to conjecture that such a property holds true, 
in soluble models (in the sense discussed here) for an arbitrary number of species s and in arbitrary dimensions. 
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VIII. APPENDIX A : DEFINITIONS AND ABREVIATIONS 

In this appendix we introduce the definitions which are adopted throughout the paper. In (|l^) and in the following, 
we have used these notations 
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_ T^20 
= 01 " 




1 p22 
r i 01 


^0 ' 




p20 , 
^ 20 ^ 


-p21 
" -L 20 " 


i-r^o 


^0 ' 




^ 02 "T ^ 02 " 




A'' 

^0 


cl 


_ p02 

— ^ m 




, p22 
^ 00 


; 




J- 10 


1 p22 
1^ J- 10 


-cl 


. r^h 

, L/2 = 




T^12 

J- 01 


, p22 
+ ^ 01 






D\ 


_ p02 
— 20 






~ "-^0 > 




p02 , 

02 1 


pl2 
-L 02 ' 




-cl 













(141) 



Further we also use the following notations: 







pll 

^ 00 


; -^0 = 


- p22 . 

- 00 I 


Tpab - 
Eq = 


_ pl2 
- ^ 00 




















El" 




p21 
^ 00 


; F^-. 


- ^ 10 


^ 00 ' 




p22 
^ 10 


p22 
^ 00 
















pab 




pl2 
^ 10 


pl2 
^ 00 


; Fl'' 


_ p21 
— 10 


p21 
^ 00 


; F^s 


- ^ 01 


^ 00 














Fl 




T^22 
^ 01 


p22 
^ 00 


; Ff 


_ pl2 
= 01 


pl2 
^ 00 


; F^^ 


_ p21 
= ^ 01 


p21 
^ 00 


















pll 
^ 20 


^ 00 


■.Fl. 


_ p22 
- ^ 20 


p22 . 

^ 00 ' 


pab^ 


- pl2 

- 20 


pl2 
^ 00 














pba 




p21 
^ 20 


p21 
^ 00 


; F2 E 


- 02 


^ 00 I 


Fl^ 


p22 
02 


p22 
^ 00 














pab 




pl2 
02 


pl2 
^ 00 


; 


_ p21 
— 02 


p21 
^ 00 


. /^a - 

, Ui - 


- ^ 00 


, pll _ 

1 J- 11 


- 

^ 01 


^ 10 










G'l 




p22 
^ 00 


1 p22 
' -L 11 


p22 
01 


p22 
^ 10 


. /^ab 
, Lri 


_ pl2 
— -'^ 00 


, pl2 
1 -L 11 


pl2 
-L 01 


pl2 
^ 10 


/^ba — 
- 


_ p21 
- ^ 00 " 


I -L 11 


p21 
" ^ 01 " 


p21 
" ^ 10 


Ci 




pll 
^ 00 


+ ^\\ 


^ 02 


-L 20 


. /^b 
■ ^2 = 


_ p22 
- ^00 


p22 
^ 02 


p22 
^ 20 


p22 . 

^ 22 1 


Gf = 


pl2 I 
-■^ 00 + 


pl2 
^ 22 


pl2 
^ 02 


pl2 
^ 20 


G\- 




p21 
^ 00 


+ ri 


p21 
^ 02 


p21 
^ 20 


; 


- -"^ 21 


^ 00 


^ 01 


pll . 

^ 20 ' 


H^ = 


pll 1 
^ 00 + 


^ 12 


^ 10 


pll 
02 


Hi 




p22 
^ 21 


1 p22 
+ -'^ 00 


p22 
^ 01 


p22 
^ 20 


; = 


_ p22 
- ^ 00 


I p22 _ 
r J- 12 


p22 
^ 10 


p22 . 

^ 02 J 


Hf^ 


pl2 , 
00 ^ 


pl2 
' J- 12 


pl2 
" ^ 02 


pl2 
" ^ 10 


Hl'^ 




p21 
^ 00 


, p21 
+ -L 12 


p21 
-L 02 


p21 
^ ^ 10 


; ^ff 


_ pl2 
— 00 


, pl2 
+ ^ 21 


pl2 

^ -L 20 " 


pl2 
-L 01 


; Hl^ ^ 


_ p21 
- ^ 00 


, p21 
+ ^ 21 


p21 
-L 20 


p21 
-L 01 



(142) 



Notice that in the expressions above, (141) and (142), the rates F^^ have not been made explicit for brevity. 
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IX. APPENDIX B: DEFINITIONS FOR SECTION V 



In section V, in addition to (141) and (142), we will also use the following additional definitions 



Aa 
Ba 
Ca 

jyA 



jyA 

^1,1 



V. 



2,0 

A 
24 



2(Ag + C„") 
2(Al + C^) 



AAd ( Ba 



— Ca+Aa] dpAioo), 



^ -Ca + Aa] d{pA{oo) - pa{0)), 



(Ff 



Aa 



Aa 
2 



Ca)pa{^) + {F^ + F^)pb{^), 



(Fi + FS - ^-f~CA){pA{Q) - Pa{^)), 
{F^ + f:){pb{Q)~pb{^)), 



(143) 



and 



As 
Bb 
Cb 



-r,B 

^2,0 



V. 

V. 



2{Bl + Dl), 



B'2 + 

ABd 



-Cb + Ab] dps(oo), 



-Cb+ AbJ d{pAioo) ~ /9a(0)), 
E'o + (F^ + F^)pAi^) + {F^ + F^-:^- Cb)pb{^), 



- 77fc 



B _ 
2,1 

B 



Fl- 



Ai 



-CB)ipB{0) - Pb{oo)), 



^(F^ + F^){pa(Q)-Pa{^)), 



(144) 



Aab,i — Ab/2, Aab,2 — Aa/^ 
C 

jsAh 

■^1,0 

■^1,1 



AB = ^2 I -^1 ^ ^2 ' ^1' 
AS _ ( aAB „ / „N , aAB 



D\ 



Bl + D^, 



- {At"" Pa{oo) + A^'' Pb{oo)) d 

-At"" ipAiO) - PAioo)) d 



-4^ {pB{0)-pB{^))d 



jyAB 

V^^ ^ E^" + {Ff + - C'M^) + [Ff + Ff - ^g)ps(oo), 



AB _ 



,AB 
2,2 



« + -Co^)(p^(0)-p^(cx3)) 
[Ft + Ff-C^,){PB{Q)-PB{^)). 



(145) 



Where pa(oo) 



and ps(oo) = 



2|7B 



as in section III. 



X. APPENDIX C: ONE DIMENSIONAL TWO-POINT CORRELATION FUNCTION ON V2, THE 
EQUATIONS OF MOTION AND THEIR SOLUTIONS 

In one dimension the equations of motion (|9^p8|) of the correlation functions can be written as an unique difference 
equation: 



25 



+ (l?^o + Kl^''^') ^rfi + (Ko + Kl^''^' + T^2,2e''^') i^r,l + Sr,-l) 
+ [(G? - BA/2){Sr,l + 6r,-l) ~ 2CA6r,o] Gl^it) 

The solution is 

+ f dt'e-I^^K*-*') {V^^Ir{2CA{t - t')) + I?2^o [Ir+l{2CA{t - t')) + Ir-l{2CA{t - t'))]) 

Jo 

+ / di'e-l^-^l(*-*')e-l'^-^l*7^(2CA(t-i')) 
Jo 

+ f dt'e-l^-^l(*-*'H2?2''ie^'^^'*' + 2?2^e"'^"'*') [/r+i(2CA(t - t')) + /._i(2C^(i - t'))] 
Jo 

+ [ dt'e-l^^l(*-*')gi^^(i') [(G? - B^/2) (/,+i(2CA(t - t')) + /,-i(2C^(t - t'))) - 2CAlr{2CA{t - t'))] 



Similarly for the B — B correlation functions we get 

j^Qr'^it) = Bsgnt) + Cb + g^^^At)) + AsPBit) 

+ [{Gl - BB/2){6r,i + Sr,-i) - 2Cb<5.,o] gi^'it) 
and for the A — B correlation function, we have 

J^gr'^it) = BABg^^{t)+CAB {gf^iit) + g^^,{t))+AABAPA{t)+AAB^2PBit) 

+ + i7f - {Al + Z?^))((5,^i + <5.,-i) - 2CAB5rfi] g^^{t) 



(146) 



(147) 



(148) 



(149) 



Eqs. (148) and (14£) are solved in asimilar way as (147). The above expressions for gi?{t), £ {A,B) can be 

rewritten in a more compact form (|99|-101) using the properties of the modified Bessel functions In{z). 



XI. APPENDIX D: CORRELATION FUNCTIONS ON V2 IN ARBITRARY DIMENSION. THE 
EQUATIONS OF MOTION AND THEIR SOLUTIONS 



The equations of motion are the higher dimensional counterparts of the previous equations ( 146| ,148,149), i.e., 
d 



and 



^qAA 

^^y\r\ = \(ri,...,ra) 



(t) = BAdgf^it) +CaJ2 {Sri'i*) + (^)) + dAAPAit) 

a 



a=l...d 



{GI-Ba/2)Y,{5. 



a' 



a^l...d 



(150) 



I n '^'■.'^ 
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{GI-BbI2)Y,{6.^ 



ol-^ol' OL — \...d 



(151) 



and also 



^y|.| = |(n,...,..)|W = BABdg^''{t)+CABY.{Q^?(t) + Qt^{t)) +dAAB^PA{t)+dAAB,WB{t) 

OL 

a=l...d 



(152) 



a=\...d 

The solution of ( |150| ) is : 



■ / di'e-l^-^l'*(*-*')l?i^o J] /,j2CA(t-i')) 

•^0 ' a=l...d 

•^0 



[| /,^,(2CA(t-t')) 



(/,„+i(2CA(t - t')) + /r„-i(2CA(t - i'))) 



+ /'*d<'e-l^-^l'*(*-*')l?^e-l^^l''*' TT /,„(2CA(t-t')) 

+ di'e-l^^l'^(*-*')g|^|ii(t')(G? - i?A/2) V 
Jo 



/,^,(2CA(t-t')) 



[| /,^,(2CA(t-0) 



(/.„+l(... )+/.„-!(...))) 
(/,„ + !(...) +/.„-!(...))) 



2CAd ( di!e-\''^\''^'-'"^Qtt^(i!) TT e-l^-^l'^(*-*')/,j2CA(i-i')) 



(153) 



where the abreviated notation (. . .) instead of (2CA(i — i')) been used. Other correlation functions QyJ^it) and 
5j^j^(t) are obtained in a similar way. Properties of the Bessel functions and elementary manipulations lead to the 
more compact forms (121-123). 
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